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I. INTRODUCTION 



Inflation [1| is known to be the best solution to the fine tuning problems associated with 
the initial conditions of the universe in the standard big bang cosmology: it resolves so-called 
the "homogeneity problem" and "flatness problem." Moreover, quantum fluctuation of the 
inflaton flelds generated during inflation can eventually provide seeds of large scale structure 
(LSS) observed in our universe Thus, it has been tested by observing the fluctuations in 
the LSS, the temperature anisotropy in the cosmic microwave background (CMB), and so on 



through the satellite experiments such as COBE [3| and Wilkinson Microwave Anisotropy 
Probe (WMAP) so far, and also will be more precisely checked by Planck in the near 
future. 

The recent WMAP seven data 0] shows that the initial density perturbation is almost 
scale invariant and Gaussian with = 0.96l;g;Q]^3 and —10 < /nl*^' < 74 at the 95% confl- 
dence level. The single fleld chaotic inflation is well consistent with the current observation, 
implymg very small uou-Gaussian,ty suppressed by the slow-roll par ameters g If sizable 
non-Gaussianity will be detected in Planck satellite experiment in the near future, however, 
it will play the role of the criterion to select a realistic inflationary model. 

Inflationary scenario is based on scalar fleld theory. Although inflationary scenario seems 
to be inevitable in cosmology, it is very non-trivial to realize in quantum fleld theory. In order 
to keep the small inflaton mass against quantum corrections, introduction of supersymmetry 
(SUSY) is helpful as in particle physics. Unless an inflationary model is not very elaborately 
constructed, however, introduction of SUSY is not enough: the positive vacuum energy in 
supergravity (SUGRA) induces the Hubble scale inflaton mass, yielding r] ~ 0{1) during 
inflation, which destroys the slow-roll condition. It is called the "eta-problem." 

"Hybrid inflation" [7| is basically an two fleld model with one as inflaton and the other, 
called waterfall fleld, to terminate the inflation when it becomes tachyonic. The advantage 
of it is that the inflaton's fleld value is small compared to the Planck scale, and thus it is 
egitimate to use the low energy effective theory. In the SUSY version of hybrid inflation 
8, 0|, the potential can be made flat enough, avoiding the eta-problem: fortunately the 
Hubble induced mass term is accidentally cancelled out with the minimal Kahler potential 
and the Polonyi type superpotential during inflation. The speciflc form of the superpotential 
can be guaranteed by the introduced U(1)r symmetry. 

By the logarithmic quantum correction to the scalar potential, the inflaton can be drawn 
to the true minimum, leading to reheating of the universe by the waterfall flelds. Moreover, 
thanks to such a logarithmic correction, the vacuum expectation values (VEVs) of the 
waterfall flelds can be determined with the CMB anisotropy j^. The VEVs turn out to 
be tantalizingly close to the scale of the grand unifled theory (GUT). Accordingly, the 



waterfall flelds can be regarded as GUT breaking Higgs in this class of models 10 



13|. This 
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inflationary model predicts a red-tilted power spectrum [9| around 

n^^l + 2r]^l- — ^ 0.98 (1) 

Nf 

for Nf = 50 — 60 e-folds. It is too large compared to the present bound on the spectral 
index. 

In the "natural inflation" model [l^, the inflaton field is regarded as a pseudo Nambu- 
Goldstone boson. Hence, a U(l) Peccei-Quinn symmetry, U(1)pq, should be assumed to 
be there, and the inflaton's small mass can be protected against quantum corrections. Its 
scalar potential, which is given by a sinusoidal functional of the inflaton, can be induced by 
instanton effects, which break U(1)pq. Since intanton effects still respect a shift symmetry, 
the axion does not appear in the Kahler potential in the SUSY version of the natural inflation 



model 15|. As a result, the unwanted Hubble scale inflaton mass term is not induced in the 
SUGRA potential. 

However, for the slow-roll parameter "77" to be small enough in this model, the Peccei- 
Quinn breaking scale or the axion decay constant / must be larger than the Planck scale, 

/ > 3Mp , (2) 

where Mp is the reduced Planck mass (~ 2.4 x 10^^ GeV). It implies that U(1)pq should 
be valid above the Planck scale. However, such U(1)pq is not natural, because quantum 
gravity effects are known to break all global symmetries including U(1)pq.^ 

In this Letter, we attempt to improve such shortcomings in the two inflation models, 
SUSY hybrid and natural inflation models, by combining them. That is to say, we will 
examine the possibility to achieve the spectral index of 0.96 and / <^ Mp in the combined 
model. We will call it "natural hybrid model" for inflation. Accordingly, we have two 
inflaton fields in this model in addition to the waterfall fields. 

One of the important properties of two or multiple fields inflation models is the existence 
of the non-adiabatic perturbations during inflation and they can change the evolution of 



the curvature perturbation after horizon exit 18l-l21|. The residual isocurvature perturba- 



tion may cause an observably large non-Gaussianity |22|, and becomes adiabatic. After 
the isocurvature modes are exhausted, the curvature perturbation is preserved during the 
radiation era and finally leaves the fossils in the CMB anisotropy observations. 

Even though it was in the slow-roll phase, the non-adiabatic mode can redistribute the 
shape of the perturbation during inflation and finally leads to the non-Gaussian feature at 
the end of inflation. The general condition for the large non-Gaussianity was analytically 



derived in |23|, 
inflation model in 25 1 



24l | with some examples and studied in detail for the multi-field hybrid 



^ It might be possible to obtain the effectively large / from sub-Planckian Peccei-Quinn scale though with 
multiple axion fields 0, [l^l ■ 
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There are other multiple field models which generate large non-Gaussianity such as the 
curvaton scenario, modulated (p)reheating, and an inhomogeneous end of infiation. Inter- 
estingly enough, it can be shown that these models have common features in the mechanism 
generating large non-Gaussianity J26| . The on-going projects to improve the observational 
sensitivity with Planck satellite |5| and using LSS data are expected to greatly constrain 
the non-Gaussianity parameters or discover the non-Gaussian feature in the primordial den- 
sity perturbation in the near future. They could discriminate a right one among various 
mechanisms of density perturbation generation. 

This Letter is organized as follows. In section II, we will construct an inflationary model 
and check the SUGRA corrections of the model. In section III, we will discuss the various 
features, particularly the spectral index and large non-Gaussianity predicted in this model. 
We will conclude in section IV. Finally, section V is devoted to Appendix, in which we 
presents various expressions related to the spectral index and the non-linear parameter, 
based on 6N formalism. 



II. THE SUGRA MODEL 

To preserve the small inflaton masses during inflation against the Hubble scale SUGRA 
corrections, let us introduce the U{1)r and a shift symmetry. Under the U(l)j^, symmetry, 
the superpotential W and a superfleld S are supposed to undergo the same transformation, 
i.e. W — > e^^^W and S — )■ e^^'^S. Under the shift symmetry, a superfleld T is supposed 
to transform as T — )■ T + 27rz/, where / is a constant with mass dimension one. We also 
consider the superflelds of a conjugate pair, tp and tp, which are assumed to carry opposite 
gauge charges. 

The superpotential consistent with the U(1)r and the shift symmetries is written as 

W = kS [M^ -m'^e-^^^ + pe-^/f)] , (3) 

where M^, m?, /, k, and p are parameters. Actually, e'^^^-^ , e^^"^/-^, etc. also can contribute 
to the superpotential Eq. ([3]). However, they do not give rise to qualitatively different 
patterns of inflation, compared to a simple case only with e~'^^^ . So we will neglect them for 
simplicity. The "p term" (p <^ 1) does not affect the inflationary scenario, since (ipip) = 
during inflation, as will be discussed later. On the other hand, it can be important when 
(ipip) 7^ 0. We assume a hierarchy among the dimensionful parameters, m f M 
Mp). It turns out to be necessary for slow-roll of the inflatons and large non-Gaussianity. 
For simplicity, we also assume that they all are real parameters. 

We note that the U(1)r forbids S*^, 5''^, etc., which destroy slow-roll inflation, and restricts 
the superpotential to the linear form of S. This superpotential keeps the small inflaton mass 
against the large SUGRA correction of the Hubble scale during inflation, as will be seen 
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later. The scalar component of T is composed of two real scalar fields, T{x) = 4>{x) +ia{x).^ 
Because of the shift symmetry, the Kahelr potential, which is a functional of T + T*, does 
not contain a{x). Accordingly, the F-term scalar potential in SUGRA is expected not to 



induce the Hubble scale mass term for a{x) during inflation 15| 
From Eq. ([3]), we can obtain the F-term scalar potential: 



Vf/^^ = \M^ - m'e'^l^ - ^^,{1 + pe~^l^f (4) 



— ,2 6 



-24>/f 



P 



The SUSY minimum is located at S = a = and — m^e"^/-^ — + pe~^^^) = 0. 
The D-term potential constrains the VEVs of and to satisfy = We assume 
that they develop VEVs in the real direction: {ip) = {^pr) / = {ip) = (^r)/v^5 where the 
subscript "r" indicates the real component for each field. At the minimum of Eq. (j4]), then, 
the VEV of a vanishes. Since we have only the two constraints in the field space {ip,ip,(p), 
i.e. M2 - m^e-'^/f - V^r?r-(1 + pe"^'f)/2 = and ^pr = '^r^ the VEVs of and ^r'^^ are 
not fixed yet, and so there potentially remains a modulus. However, by including the soft 
SUSY breaking mass terms for them, dVsoit = {'m'^(p'^ + mlipl + m'^ip^)/2^ their VEVs can 
be determined, 

U''^^-2'^'%^0{mVP) « 1, (5) 
^ 2{M^ - m^e"^'^ - ^-^mj) ^ 2{M'^ - - ^-^mj) , (6) 

where m^, m^, and m:^ indicate the soft mass parameters for the real components of </>, ip, 
and ip. They are of the TeV scale. We set = m-^ for simplicity. From Eq. ([5]), we see 
(p/f^ — 2(m^m^)/(m^/^) ~ 0{m'^ / p) ^ 1. Around the minimum of the scalar potential, 
the mass eigenstates {$, \E'} and their corresponding mass squared turn out to be 

(p + e^r , ml^ml ; and ^ ipr - e(p , m| ^ 4k^M^ , (7) 

where e = m^e-'*'/ ^ / y/2f M . Thus, the VEVs of $ and ^ are 

($) ^{(p + eVv) ~ O [m^/f) , and (^) ^ {ip., - ecp) ^ V2M . (8) 

We will focus on the case m^// ^ 10^° GeV, because this case turns out to yield large 
non-Gaussianity. 

To obtain a sufficient inflation, we suppose S > M initially. Since jS'p plays the role 
of the mass squareds for ip and ip, the initial condition S > M compels the VEVs of ^p 
and ip to vanish during inflation. As a result, S and T can be light enough and the scalar 



In this Letter we use the same notation for a superfield and its scalar component. 
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potential becomes dominated by the positive vacuum energy density k^M'*, only if the 
mass of S in Eq. is much lighter than the Hubble scale, i.e. /t^m^//^ <^ n'^M^/Mp or 
m^/M^ <^ /"^/Mp 1). Thus, the condition S > M provides a quasi-flat scalar potential 
with positive vacuum energy density. Under the condition, thus, inflation can arise. 

The positive vacuum energy density during inflation breaks SUSY explicitly. The coupling 
between S and the waterfall fields in Eq. ([3]) and mass splittings between the scalars 
and fermions by SUSY breaking induce the quantum correction at one loop to the scalar 
potential Eq. (Q: 

SVf ^ X ^ logf , (9) 

where /i^ (= k^M^) denotes the vacuum energy density during inflation and A the renormal- 
ization scale It can draw the inflaton S to the true minimum of S* = 0. Even with this 
logarithmic functional of S, the slow roll conditions still holds. We assume it is dominant 
over the mass term of 5* proportional to / p in Eq. (jl]). 

As 5* approaches to M , if) and if) become destabilized, since they become tachyonic for 
S < M around the origin. Thus, ip and ip start rolling down to their true minima discussed 
in Eq. (|4]), when S becomes smaller than M. It means that the VEVs {ip) and {ip), which 
play the role of the masses for S as seen in Eq. (jl]), grow from zero approximately up to 
M. We note that before ip and ip reach M, the slow-roll condition for S breaks down when 
r]s = MldlV/V ^ 1, that is to say. 

Ml + m')) ^ V , (10) 

where V is approximately given by Eq. (j4j). This is the condition for end of inflation in this 
model. Since the leading term of is a large constant (k^M^), the left-hand side of Eq. (fTOjl 
should be so. For instance, {ip) (= {ip)) is approximately given by M'^ /Mp !^ 10"'^'^ GeV, if 
M ^5x 10^5 GeV. 

With the D-flat condition, \^\ = \^p\ = \^r\/V^, Eq- ( ITOll is written in terms of 
MIk'^^PI ^ K^[M^-{M^-\S\'^)^Pl + %, where ^ M^(l -2|^cos|), and we neglect the 
term proportional to m^. Since the waterfall field ipr is sitting always at the local minimum, 
its VEV is determined by d^^V = 0, yielding ^ = - \S\^ = ^/M^ - V/k^. It recasts 
Eq. (do]) into +Ak^M^V -Ak'^MI,M^ = or 1/ = 2K^Mj,{-l + y/l + M^Mj,) ^ k^M^, 
i.e. V{S, a) is almost a constant (~ k^M^) when inflation is over. It implies that the 
condition for end of inflation almost respects a uniform energy density condition in this 
case. Hence, one can expect that the end point effect in non-Gaussianity would be quite 
small Uj 



In Rcf. 



281 ] ■ the end point effects in non-Gaussianity, particularly for the cases that two inflaton masses 
are degenerate (mf = to|) and hierarchical {ml ^ m|), have been studied. Following the notation of the 
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Once the slow-roll condition for S violated, that of a is also violated quickly. During 
— tjjip ^ + pipip in Eq. (jl]), the slow-roll parameter of axion T]a (= Mpd^V/V) can be 
smaller than the unity even with / -C Mp. As — ipip decreases and becomes comparable 
with or pipip, however, the slow-roll conditions for a becomes also violated, since rja 
becomes Mp/ f'^ 1) in this case. Thus, a should also roll down to zero. Due to the 
existence of the non-adiabatic mode, the total curvature perturbation evolves until the end 
of inflation and abruptly stops to evolve after slow-roll is violated. The perturbation of 
waterfall fields during the waterfall dynamics does not change the curvature perturbation in 
the large scales as studied recently in many literatures |30l-l35[|. 

Finally, let us check SUGRA corrections. We ignore the quantum correction Eq. <^ for a 
while. By including SUGRA corrections, the Hubble scale mass terms for some scalar fields 
could be induced during inflation. To see which fields acquire such masses, let us consider 
the full F-term scalar potential in SUGRA Lagrangian. Since if) and if) got already heavy 
masses proportional to S even in global SUSY, we will not consider them. 

During inflation, thus, the superpotential and the Kahler potential are given by 

W = kS{M^ -m''e-^/f) , K = \S\'' + hj + T*f , (11) 

respectively. They respect the U(l)/j and shift symmetries. Here, the Kahler potential takes 
the minimal form without containing T — T*. The resulting kinetic terms of the scalar fields 
are of the canonical type. We will see no Hubble induced mass term for S appears in the 



scalar potential (36|]. If there exists a quartic term of S in the Kahler potential, however, 
it would generate the unwanted Hubble induced mass term, which destroys the slow-roll 
conditions for S. Hence, we should assume its coefficient is small enough (< 10"^), if it 
exits. This assumption needs to be justified by a proper UV theory of SUGRA for its 
naturalness in the future. Since the shift symmetry is encoded in the Kahler potential, only 
the real part of T, i.e. (f){x) [= [T + T*)/2] appears in it. Due to the reason, we will see 
that the Hubble induced mass term for a{x) [= {T — T*)/2i] does not appear in the SUGRA 
potential. 

With the covariant derivatives in SUGRA, 

^ dW W dK 2 A \S?\f-, r/A 

DrW^?^ + iil^ « n-I^ (t + Ml^-^eM (13) 
^ &T Ml&T Ml \ J P KP ) • ^ ' 

Ref~j3|, G = 2k^|S'P and = 2k^A4^, respectively, in our case. Thus, our model corresponds to the 
limit, gi=a = 6~0, 7r/2 of Ref. [28|, and mf and m| are estimated as ria\s=M and rjsls^M, 

respectively, and as will be seen later, rya|s=-A/ S> ?7s|s=a/ in our case. With such limits of the parameters, 
the condition for end of inflation corresponds to a straight line in the field space, and the end point effect 



in non-Gaussianity can not be large as noticed in 27|,|29|. It is consistent with the analysis of Ref. [28 



7 



one can write down the F-term scalar potential: 



UGRA — e ' P 



9 9 \W 

\DsW\^ + 2\DtW\^ 



Ml J 



-T/f 



Ml 



Ml 

7 ^ /2 M2* 



(14) 



where we used e^/^'^p ^ (1 + 0VM|,)(1 + |^|VM|,), and i^^^* = 1/Kss* = 1, /^^ 
1/KxT* = 2, -R'sT* = 0, etc. We dropped \S\*/Ml in Eq. f lT^ . because of its smallness. 
Due to the presence of the last term, fi^(f)'^/Ml 3if^</)^) in the last line of Eq. f|T^ . which 
spoils the slow-roll condition for (f), (p should be stabilized somewhere during inflation. On 
the other hand, the quadratic term of S was cancelled out, and so it does not appear in 
Eq. ([Tlj). It is a nice feature of the SUSY hybrid inflation model Since "e^/^p" in 
the scalar potential does not contains a, the Hubble induced mass term of it also does not 
appear in Eq. (fT4|) . It results from the shift symmetry 151]. 

In our study we require the hierarchy between the parameters 



m 



Ml 



< 1 



(15) 



We will see later that this condition is necessary for slow-roll of a. Note that / is re- 
garded as being smaller than Mp. In this case, the VEV of (p is estimated as \4>/f\ ~ 
0{Mlm^/pM'^) < 1. So we can neglect "e-<^/-^" and "0V/^" in Eq. ([14]). 



III. NATURAL HYBRID INFLATION 

We are particularly interested in the following form of the potential for inflation with two 
scalar fields: 

Mnf = {l + alog^ - -^cosy^ , (16) 

where the scalar fields S{x) and a{x) play the role of the inflatons. It can be derived from 
Eq. (jlj) or (fT4|) by taking ?/; = '?/' = 0, 0//^l, and Eq. (fT5|) . The dominant vacuum energy 
comes from the first term /i^. The second term in Eq. (fT6|) is originated from the scenario of 
the SUSY hybrid inflation and third term from the natural inflation. As discussed in Eq. ([2]), 
the second term is generated by SUSY breaking effects at one loop ^ and so a = k^/Stt^. 
We assume a hierarchy between the dimensionless parameters a and A, 1 ^ a ^ A. Thus, 
the dynamics of the trajectory in the field space is mostly dominated by S field, though 
the cosmological observables such as the spectral index and non-linear parameter /nl are 
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controlled by both fields. By comparing Eq. (fT6|l with Eq. (fT^ . the parameters /i'* and A 
can be identified as 

,^ = ^M\ and A.^e-^//-^. (17) 

With Eq. ( IT6l) the slow-roll parameters for S and a are estimated as 

_ a^Ml _ a _ aMl _ 1 

where the fields and parameters were simplified as 

Hence, the slow-roll conditions (e^^a <^ 1, |?7s,a| ^ 1) are fulfilled as long as a and A are 
small enough and 

» 1 , and = ^A « 1 , (20) 

where we see that the second condition is consistent with our previous requirement Eq. (fT5l) . 
Note that if A is quite small, "/" can be much smaller than Mp unlike in the original natural 
infiation scenario. 

During the slow roll infiation the equations of motion for 5 and a fields are simplified to 

be 

4 

3HS + ^ + An^ij^S = 
3if^ + ^sin^ = 0. 

Neglecting the ip term, which vanishes during infiation, they give the solutions for the 5* (or 
X) and e (= a//), 

Xl - xl = 2Nf , and tan| = tan| e^'^^ , (22) 

where Nf is the e-folding number between the horizon exit and the end of infiation, which 
is assumed to be around Nf ^ 50 — 60. Since we require slow-roll condition until the end of 
infiation, we take xl > 2iVj ~ 100. Throughout this Letter, the subscripts (or superscripts) 
"*" and "e" denote the values evaluated at a few Hubble times after horizon exit and the 
end of infiation, respectively. Since infiation is over when S'e ~ M (= /i/^^), we have the 
relation, V\e ~ /i^ ~ Sn'^a^Mpxt- 
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Since we are interested in the evolution of the curvature perturbation in the superhorizon 



scale, we use 6N formalism 



37-41 



The number of e-foldings, A^, is given by 



N 



H{t)dt, 



(23) 



and the curvature perturbation is evaluated as 



(24) 



IJ 



where N,i = ON/ d(fl and the index I runs over all of the inflaton fields. We will consider 
the power spectrum and bispectrum defined (in the momentum space) by 



27r2 



(25) 
(26) 



With the 6N formalism the power spectrum, the spectral index and the non-linear parameter 
are given by 



4^ 



— 1 



-2e* 



2 Eij^jNjiNi 



K 



fl 



NL 



6 {Y.iN-jf 



(27) 
(28) 
(29) 



if the infiatons' kinetic terms are given by the canonical form. Especially for the separable 
potential we can calculate e-folding number analytically and so the analytic expression of 
the corresponding /nl can be obtained 



23 



42| . For sum potentials, we summarize them 



in the Appendix. Therefore, the curvature perturbation evolves after horizon exit during 
the inflation and the non-Gaussianity can be developed, which can be observed through the 



satellite experiments in the near future 2J, |25|. More examples of large non-Gaussianity 



43 



45 



471. 



generated during the inflation with the adiabatic limit were studied recently 

Field fluctuations are generated at the horizon exit with almost Gaussian statistics, and 
thus the curvature perturbation is also Gaussian at that time. However, due to the non- 
adiabatic modes in the multiple field inflation, the curvature perturbation experiences non- 
linear evolution and changes its shape to be non-Gaussian. When the non-Gaussianity 
becomes larger than 0{1) during inflation, we can approximate it analytically as in the 



papers [2J, |25[. The size of the non-linear parameter /nl is determined by the slow- roll 
parameters e^'s and r/j's. Especially the sign of /nl is given by the relative size of rj between 
the horizon exit and the end of inflation by —77* + 2?7e of a relevant field 24 1. 
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In our model, we have two inflaton fields S and a, which are slowly running during 
inflation with the waterfall flelds, ip and ip flxed at the origin. Large non-Gaussianity can 
arise when the axion fleld a is located around the top at the horizon exit. Since the axion 
fleld is on the ridge where the potential is concave, i.e. r]a < 0, the trajectories diverge and 
large negative non-Gaussianity is generated. After the axion fleld crosses the convex point, 
a/f = 7r/2, the potential changes the curvature and rja > 0, which makes the change of 
the sign of /nl-^ Soon the waterfall fleld is destabilized and the inflation get to the end, 
when the slow-roll condition is violated. However, the waterfall dynamics does not affect the 



curvature perturbation after slow- roll inflation as studied in Refs. |30l-l35j . The dynamics of 
the axion fleld and the evolution of /nl are shown in Figs. H] and [2] respectively. Here we 
used the parameters of Case 1 in Table [11 

In the region of large non-Gaussianity, we can make good estimates for the power spec- 
trum and non-Gaussianity using the analytic study in the separable potential in Refs. 24|.|25|. 
In this region of large |/nl| with ^ ea, we get 

SlN^e = ^(-] sin^O « 1 , (30) 

+ ea V"/ 

where we used Eq. ([T8|) . Therefore the power spectrum (P^), spectral index (n^), and 
tensor-to-scalar ratio (r) are approximately given in terms of the slow-roll parameters: 

^^-24^^^ + ^-)' (^1) 
nc-l--2e: + 2Z^l±|^, (32) 

(33) 



1 + f ' 



where f is deflned as the ratio of the contribution to the curvature perturbation from each 
fleld by 



Ns SIN'e, " \aj xl sin'e. 

In Appendix, the relevant original expressions are summarized. 

In the limit of f — > 0, the linear term of the curvature perturbation is dominated by 
S fleld and the power spectrum and the spectral index is determined solely by the S fleld 
perturbation. Thus it reproduces the result of the SUSY hybrid inflation, i.e. ^ 1 + 2//* ^ 
0.98 from Eq. (ITSl) and Eq. (l32l) with xl > 2A^/. In the opposite limit, for large f ^ 1, the 
power spectrum comes dominantly from the a fleld perturbation and, the spectral index is 

nc - 1^27/: = 2^2 cos ^ , (35) 



^ The behaviors of /nl on the ridges and valleys are also well described in the paper of EUiston et al. [45 1. 
See also Ref. 
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FIG. 1: Evolution of axion field without waterfall field (with blue line). However, the inflation 
ends at e-folding number Nf = 60 (with vertical dotted line) and the axion decouples from the 
blue line and moves quickly to the minimum of the potential. 

f_NL 



200 - 



100 ^ 



-100 L 




FIG. 2: Evolution of /nl for case of a positive /nl (Case 1 in Table H]). The dashed line corresponds 
to the case without waterfall fields, and the solid line is for the model with the waterfall fields, 
where end of inflation occurs at Nf = 60 (with dotted line). 

so that we can chose ^'^ ~ 0.02 to make ra^ ^ 0.96 for a^,// ~ tt. To get large positive 
/nl, however, we need > 0.1, because we need 9e < vr/3 from /nl oc —t]* + 2r]^ oc 
— cos 6^ + 2 cos ~ — 1 + 2 cos 6e, and equation of motion Eq. ( I22l) 



Oe 2 e-«'^/ 

tan ^ ^ — . (36) 

2 TV — 

This predicts 

nc < 0.8 , (37) 
which is not consistent with WMAP results. Therefore, the only available value for f com- 
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Case 


^^{= 




a (= 




A(« 


2mVM2) 


/ 




TT 


— 0^ 


e 


X* 


r 


/nl 


1 


1.7 


1014 


2.6 


• io-« 


5.0 


• 10-13 


5.4- 


1012 


1.3 


10-2 


0.10 


15.9 


0.19 


43 


2 


1.7 


1013 


8.6 


10-11 


2.4 


• 10-17 


3.8 • 


lOio 


1.5 


10-2 


0.10 


86.1 


0.24 


72 


3 


9.3 


lOii 


7.8 


10-14 


2.0 


• 10-22 


1.1 


10^ 


1.5 


10-2 


0.096 


887 


0.23 


67 



TABLE I: Some parameter values which give large /nl at the end of inflation with Nf = 60. The 
values of fi and / are listed in the unit of GeV. Here we imposed the constraints from Power 
spectrum Eq. (I38p and the spectral index Eq. (I39p . M (= [ij ^/k) and / should be around the GUT 
and intermediate scales, respectively. 



patible with large positive /nl and the spectral index is f = For small /nl, however, 

there should be no problem even with much smaller ^2. w^e listed some of cases which can 
predict large non-Gaussianity in Table [H In Fig [3], we show the evolution of /nl for a case 
with negative non-Gaussianity. 

In Table [H we imposed the data from seven years of WMAP 

= 2.43 ±0.115 X 10-\ (38) 
= 0.96lJ^:[J}3 (assuming r = 0) . (39) 

By comparing Eqs. (13T|) and (!32|) with Eqs. (!38|) and (!39|) . we can get the relations 

^a2^4^2(i^~)^ 2.43x10-9, (40) 

1^(3 + f) ^ -1 + fx'e'cos^* + 0.02(1 + r)xl • (41) 

Hence, a and f can be determined in terms of x* (or Xe), , and ^2 Since a is a quite 
small positive number, neglecting the left-hand side of Eq. ( 14T|) . f would be 

. - VX^ + 0.02 ^^2) 



^2 - 0.02 ' 

for X* 3> 1 and 9^ ^ n. Once f is determined, (A/a) is also done by Eq. flM|) . We note that 
ris ^ 0.96 is possible in the presence of the two inflaton fields S and a, unlike in the original 
SUSY hybrid inflation model. 

The tensor-to-scalar ratio is suppressed by the slow-roll parameter e* which magnitude is 
smaller than r]* by a factor a. Therefore in this model the tensor-to-scalar ratio is negligible. 

If the non-Gaussianity is large, namely for SIN26* <^ SIN29e <^ 1, it is approximately 
estimated as 24 1 

^ ^QTM^n'^n ^-^2 ("^^ + ^ afM ^ 2^-1. n ^ -^2 ("^o^^* + 2cos^e) • (43) 
6 SIN Be (1 + r)^ 6(A/a)Xe sm 6^e(l + r)-^ 
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FIG. 3: Evolution of /nl for case of a negative /nl- The dashed hne is without waterfall fields and 
the solid line is for the model with the waterfall fields, where end of inflation occurs at Nf = 60 
(with dotted line). 

Since Xe, ^e, A/a, a and f can be given in terms of x*, and by Eqs. (l22l) . (134|) . (140|1 
and (HTi) . /nl can be determined only if x*, 6**, and are given. ^ The condition for large 



non-Gaussianity |/nl| ^ 1 reads [24 1 



1 I > 1 , where Q, = ^ {-vl + 2^^' ^ . , i • (44) 



SIN'Ge^ / ' 5 I /a. /a I 5 | - COS^, + 2cOS^e 



From this condition, the following three necessary conditions can be derived 



24 



" a ) I — cos6'* + 2cos6'e | ' ^ ^ 

SIN^e, °' x^sin^^, ei - cos^. + 2cos^e| ' ^ ^ 

They are useful to search for the conditions for large non-Gaussianity. We see that large 
non-Gaussianity is possible, only if \/a and n — 9^ are small enough. For 9^ ~ vr, thus, a 
positive large /nl requires 9e < 2ti /?>. However, we don't discuss here how the inflaton a is 
initially set on the top of the potential (6** ^ vr), since it is beyond the scope of this Letter. 
For the maximum and minimum values of /nl in the 95% confidence level of the WMAP 
7-years data 0], 

- 10 < f^^l""^ < 74 (95% CL) . (48) 



^ Once Xe, a, A are known, then the SUGRA parameters k, M^, m? , and / can be determined via a = 
kVStt^, xl = Sl/aMl « AP/aMj,, X « 2m^/M^, and = M^X/f from Eqs. 1^, and (HH). 
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FIG. 4: Variation of /nl for different initial positions of the axion field with the same parameters 
used in Case 1. 



When the axion-Uke inflaton field a is around the top of the potential, i.e. 6 ^ it, 
quantum fluctuation of the field would become important and much affects its motion. For 
dominance of the classical motion, the required condition on the background trajectory is 



a\'K /H"^ > \/S/2 25|, |48]. This condition constrains the initial value of a (= fO^). Using 



Eq. (1211) . it can be recast as 

sing, > ^ ■ (49) 

The initial values listed in Table |T] fulfill the condition. 

For large /nl, M (= ji/ y/R) and / should be around the GUT and intermediate scales as 
seen in Table [11 As discussed earlier, M is the VEV scale oi ip, ip after inflation over, and 
/ could be the PQ symmetry breaking scale. Thus, the reasonable range of the parameters 
can yield the values required in particle physics. 

In Fig. m we show /nl at the end of inflation with different initial values of the axion. 
The large positive /nl is possible for the axion located around the top, (6'^, < n). However, 
it becomes negative if the axion is so close to the top. 



IV. CONCLUSION 

In this Letter, we have constructed an inflationary model by combining the conventional 
SUSY hybrid and natural inflation models to achieve ~ 0.96 and / <^ Mp. ^{1)r and 
a shift symmetry together with the minimal Kahler potential are essential to avoid the eta- 
problem. However, the SUSY hybrid inflaton sector needs to make a dominant contribution 
to the vacuum energy density during inflation. The non-adiabatic mode from the dynamics 
of the two canonical fields during inflation are responsible for large non-Gaussianity and 
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small tensor-to-scalar ratio. The symmetry breaking scale by the waterfall fields and / in 
this model should be around the grand unification and intermediate scales, respectively, for 
large non-Gaussianity. 
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V. APPENDIX 



W{^,x) = U{^) + V{x)- Defining 



u 



with 



Vp COS 9p — Up sin 9^ 



(50) 



(51) 



the power spectrum and spectral index are given by |42|: 



- 1 = -2e* - 4 



The non-linear parameter /nl is 421 ]: 



(52) 
(53) 



/] 



NL 



1 —U 



e* 



1 -V 

ez \ 2e* 



^1 



U V 



A, 



(54) 



where we define 



f, = M.P + g^^x) =^^sin2e + r/^cos^e 



A, 



y X r e _ 

(e-)2 i I 



^ cos^ Ge sin^ Ge [e" - r)^ 



(55) 
(56) 
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The slow-roll parameters are 

Ml f d^U V _ 2^ _Ml ( 



ecos^e, e^ = -^[^^^\ =esm2e, (57) 



and 



''^^^^ZTtV' ^^x = 55x^ = 0, ^, = M|.^^. (58) 
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